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Abstract 

The concepts of P- and Po-matrices are generalized to P- and Po-tensors of even and 
odd orders via homogeneous formulae. Analog to the matrix case, our P-tensor definition 
encompasses many important classes of tensors such as the positive definite tensors, the 
nonsingular M-tensors, the nonsingular H-tensors with positive diagonal entries, the 
strictly diagonally dominant tensors with positive diagonal entries, etc. As even-order 
symmetric PSD tensors are exactly even-order symmetric Po-tensors, our definition of 
Po-tensors, to some extent, can be regarded as an extension of PSD tensors for the odd- 
order case. Along with the basic properties of P- and Po-tensors, the relationship among 
Po-tensors and other extensions of PSD tensors are then discussed for comparison. Many 
structured tensors are also shown to be P- and Po-tensors. As a theoretical application, 
the P-tensor complementarity problem is discussed and shown to possess a nonempty 
and compact solution set. 
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1 Introduction 

The P-matrix, first introduced by Fiedler and Ptak m , is an important class of special ma¬ 
trices whose determinants of all principal submatrices are positive. As an essential concept 
in matrix theory (see, e.g., p~3l IT7l 18722 129]). P-matrices contain many notable matrices 
as their subclasses, such as positive definite matrices, nonsingular M-matrices, nonsingular 
completely positive matrices, strictly diagonally dominant matrices, etc., Besides its signifi¬ 
cance in matrix analysis, P-matrix also plays an important role in both the theory and the 
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computations of linear complementarity problems (see, e.g., [2] El IE] 12] ) • Given a matrix A 
and a vector q, a linear complementarity problem LCP(q, A) is to find a vector x such that 

(1) x > 0 , Ax + q > 0 , (x, Ax + q) = 0, 

where (•, •} denotes the inner product of two vectors. The matrix A is a P-matrix if and only 
if the corresponding LCP(q, A) has a unique solution for any vector q. The P-matrix is also 
a widely-applied tool in nonlinear complementarity problems. For instance, Qi et al. [23] 
employed some properties of P-matrices to investigate the convergence of smoothing Newton 
methods for nonlinear complementarity problems. 

With an emerging interest in the assets of multilinear algebra concentrated on the higher- 
order tensors, more and more structured matrices have been extended to higher-order struc¬ 
tured tensors. Very recently, Song and Qi [3T| extended the concepts of P- and P 0 -matrices 
to P- and Po-tensors. Their applications in tensor complementarity problems (TCP) (an 
extension of linear complementarity problems as studied in [31 UK GH 123]) were discussed in 
1321321. As pointed out by Yuan and You [35], there are no odd-order P-tensors according 
to the definition of Song and Qi [31]. This deficiency is more or less due to the involved 
non-homogeneous term. In this paper, we will introduce a homogenous definition of P- and 
Po-tensors, and discuss their applications in TCP and others. 

Analog to the matrix case, our P-tensor definition encompasses such notable classes as the 
positive definite tensors, the nonsingular M-tensors, the nonsingular H-tensors with positive 
diagonal entries, the strictly diagonally dominant tensors with positive diagonal entries, etc. 
Of particular concern is that our definition of Po-tensors also extend the class of positive 
semi-definite (PSD) tensors to the odd-order case. It is well-known that PSD tensors have 
massive applications in optimal control, magnetic resonance imaging, spectral hypergraph 
theory and tensor complementarity problem, etc. However, the definition of PSD tensors 
restricts the order of the involved tensors to be even. As even-order symmetric PSD tensors 
are exactly even-order symmetric Po-tensors, our definition of Po-tensors, to some extent, can 
be regarded as a counterpart of PSD tensors for the odd-order case. 

The rest of the paper is organized as follows. In Section 2, the definitions of P- and Po- 
tensors are introduced and some basic properties are characterized. In Section 3, relationship 
among P 0 -(P-)tensors, PSD (PD) tensors, and other PSD-related tensors such as (strictly) 
copositive tensors, (strongly) completely positive tensors, and doubly nonnegative tensors are 
discussed. In Section 4, some structured P- and Po-tensors are established. Properties on 
P- and Po-tensor Cones, together with their relationship to other special tensor cones are 
analyzed in Section 5. Applications of P- and Po-tensors in TCP and others are described in 
Section 6. Concluding remarks are drawn in Section 7. 


2 Definitions and Basic Properties 

The concept of P-matrix m is originally defined by that all of its principal minors are 
positive. As an important equivalent definition of the P-matrix, we know that a matrix A is 
a P-matrix if and only if for each nonzero vector x € M" there exists some index i such that 
Xj(Ax)j > 0 (see [2] ID]). Following this equivalent definition, Song and Qi [31] generalized 
the concept of P-matrices to higher order tensors and they call an m th -order n-dimensional 
tensor A a P-tensor if and only if for each nonzero vector x £ R n there exists some index 
i such that Xi(Ax rn ~ 1 )i > 0. Nevertheless, it has been shown that there are no odd-order 
P-tensors according under the Song-Qi’s definition (see 23)- To make up this deficiency, we 
now propose a modified definition of P-tensors as follows. 
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Definition 2.1 (P- and Po-tensors). An m th -order n-dimensional tensor A is called a P- 
tensor, if for each nonzero x £ R ra there exists some index i such that 

(2) x™- 1 {Ax m ~ 1 ) i > 0. 

An m th -order n-dimensional tensor A is called a Po-tensor, if for each nonzero x £ R" there 
exists some index i such that 

(3) Xi 0 and x” 1_1 (^ix m_1 ) i > 0. 

It is apparent to see that our definition is the same as the one of Song and Qi in the 
even-order case. And the expressions m and m can be rewritten as 

Xi ^ 0 and (Ax m_1 )i > 0 (or > 0), 

respectively. Thus, the properties about the P- and Po-tensors in [211SH [33] [35j are certainly 
valid for our P- and Po-tensors of even orders. Moreover, we will shortly show that most of 
those properties also holds for our P- and Po-tensors of odd orders. 

We start with the spectral property of P- and Po-tensors. According to Qi’s work in [T5] . 
for any given m th -order n-dimensional tensor A , if there is a real scalar A £ R and a real 
vector x £ R" satisfying 

Ax™' 1 = Ax 1 " 1 ” 1 !, 

then A is called an H-eigenvalue of A and x is called a corresponding H-eigenvector. Here 
x [m-i] ._ ..., x™ -1 ) . The following theorem shows the positivity (nonnegativity) 

of the H-eigenvalues of a P-tensor (Po-tensor). 

Theorem 2.1. A P-tensor (Po-tensor) has no non-positive (negative) H-eigenvalues. 

Proof. Let A be an m th -order n-dimensional P-tensor. If A has no H-eigenvalue, then the 
result is trivial. Otherwise, we use A £ R and x £ R" to denote an H-eigenvalue of A 
and a corresponding H-eigenvector, respectively. That is, Ax m_1 = Ax^™ -1 !, which implies 
x™ _1 (Ax m_1 )i = Ax^ m ^ for all i = 1,2,..., n. By definition, we know that x/0. Thus, 
we can always find some index i such that the left-hand side is positive, which further implies 
that the H-eigenvalue A must be positive as well. 

The proof of the results about Po-tensors is similar, thus we omit it. □ 


Moreover, it can be easily verified that the property of a P- or Po-tensor is also preserved 
by all of its principal subtensors. 

Theorem 2.2. If an m th -order n-dimensional tensor A is a P- or Po-tensor, then each of 
its principal subtensor is also a P- or Po-tensor, respectively. 


Proof. Let A be a P-tensor. Without loss of generality, we can assume that its first k x k x 
• • • x k subtensor Ai-.k is not a P-tensor, i.e., there exists some nonzero y £ R fc such that 
yl" _1 {A\-ky m ~ 1 )i < 0 for all i = 1,2,.. ., n. Denote x = [y T , 0 T ] T £ R”. Then we have 


Ax" 1 ” 1 


A ,, 771 — 1 

Au/ty 

* 


Hence, the products of the entries of xl m ^ and the corresponding entries of Ax m 1 are 

yr~ 1 ('W T "~ 1 )i<0> i = 1,2,... ,k, 


„m— 1/ /i m—1 


(Ax m_i )j = 


0, 


i = k + 1 , k + 2, ..., n. 


This contradicts that A is a P-tensor. Therefore, each principal subtensor of a P-tensor is 
also a P-tensor. □ 
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From this theorem, we can directly point out that all the diagonal entries of a P-tensor 
(or Po-tensor) are positive (or nonnegative), since each diagonal entry is a smallest principal 
subtensor. This is one of the most easily checkable necessary conditions of P- and Po-tensors. 

The following equivalent definition of P-tensors describes the positive definiteness of a 
P-tensor corresponding to a fixed vector. 

Theorem 2.3. An m th -order n-dimensional tensor A is a P-tensor if and only if for each 
nonzero vector x £ K™, there exists a positive diagonal matrix D x such that 

> 0 . 

Proof. For a tensor A and a vector x, denote 

/ 1, x™-\Ax m ~ 1 )<>0, 

1 \ 0 , < 0 . 

By the definition, if A is a P-tensor, then there must be some index i such that Si > 0. Thus 
it holds that E'Li (Ax m ~ 1 )i > 0. Let 

n r . r EILi SixT-\Ax^)i 

Hence, if we denote D x = diag(<5i + e, 62 + e, ■ ■ ■, S n + e), then the diagonal entries of D x are 
positive and 

n 

(D x x^- 1 \Ax m - 1 ) = Y / (^ + ^- 1 (Ax m - 1 ) i > 0 . 

2=1 

The reverse part is obvious, thus we omit its proof. □ 

3 Comparison with Extensions of PSD Tensors 

Positive semidefinite (PSD) tensors were first introduced and systemically investigated by Qi 
m , which have wide applications in optimal control (see ni). magnetic resonance imaging 
[35] , spectral hypergraph theory mmm, tensor complementarity problem Him S3 [32], 
etc. PSD tensors are a generalization of PSD matrices and are defined as follows. 

Definition 3.1 (Positive definite [Tj5][5]). An m th -order n-dimensional tensor A is called pos¬ 
itive definite (PD) if Ax m > 0 for any nonzero x £ R n ; And it is called positive semidefinite 
(PSD) if Ax m > 0 for any x £ R”. 

Actually, we can see from the definition that a PSD tensor is not necessarily symmetric. 
However, we are still most interested in the symmetric case, otherwise we can study the 
symmetrization (see [5] for details) of a nonsymmetric PSD tensor instead. Even without 
the symmetry, it is easily concluded that a PD or PSD tensor is also a P- or Po-tensor, 
respectively. 

Proposition 3.1. A positive definite tensor is a P-tensor, and a positive semidefinite tensor 
is a Po-tensor. 

Proof. If A is a zero tensor, then this is true obviously. Let A be an m th -order n-dimensional 
nonzero positive semidefinite tensor. Thus m is even. Assume on the contrary that A is not 
a Po-tensor. Then we can find some nonzero x £ K n such that 

a;™ _1 (Ax m ^ 1 ) i <0, Vi £ supp(x) := {i £ [n] : x t ^ 0}. 
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Since m is even, the above equation is equivalent to 

<0, V* £ supp(x), 

which further implies that Ax m = Y^i=i x i(Ax m ~ 1 )i < 0. This contradicts to the positive 
semidefiniteness of A. Thus A is a Po-tensor. Similarly, we can prove the PD case. □ 

The indeed incompleteness in the definition of PSD tensors is the absence of odd-order 
PSD tensors, which can be observed apparently by replacing x by —x when the order m is 
odd. Nevertheless, the odd-order cases are not allowed to be ignored. For instance, there is no 
intrinsical gap between the definitions and properties of the even-order uniform hypergraphs 
and the odd-order ones. And it is well-known that the Laplacian tensor of an even-order 
uniform hypergraph is PSD, but we have not got a parallel result for the Laplacian tensor of 
an odd-order uniform hypergraph yet. 

Song and Qi [31] pointed out that an even-order symmetric tensor is PD or PSD if and 
only if it is a P- or Po-tensor, respectively, since a P- or Po-tensor has no nonpositive or 
negative H-eigenvalues, respectively. Thereby, the equivalence between the PD tensors and 
the even-order symmetric P-tensors inspires that the odd-order symmetric P-tensors may be 
a good extension of the positive definiteness to the odd-order case. Apart from symmetric P- 
and Po-tensors, we have so far other three kinds of extensions of the PD and PSD tensor to the 
odd-order case, i.e., (strictly) copositive tensors, (strongly) completely positive tensors, and 
(strongly) doubly nonnegative tensors. We shall compare our extension with those approaches 
in this section. 

Copositivity might be the first attempt to generalize the concept of positive definiteness 
(see [SUiliEI]). Studying on copositive tensors is restricted on the nonnegative orthant M™ := 
{x £ I" : Xi > 0, i = 1,2,..., n}, which thus avoids the ill-posedness of the PSD tensors of 
odd orders and admits the PSD tensor cone as a subset. 

Definition 3.2 (Copositive [2Q])- A n m th -order n-dimensional tensor A is called copositive 
if Ax™ > 0 for any x £ MIL; And it is called strictly copositive if Ax' n > 0 for any nonzero 

x £ R+. 

For even-order tensors, the definitions directly imply that a PSD tensor is copositive. 
Moreover, we shall see shortly that this is also true for odd-order P-tensors. We need a 
lemma about the spectra of a copositive tensor from Song and Qi [30] Theorems 4.1 and 4.2]. 

Lemma 3.2 (}30]1. Let A be an m th -order n-dimensional symmetric tensor. Then A is 
copositive if and only if every principal subtensor has no negative H-eigenvalue with a positive 
H-eigenvector; And A is strictly copositive if and only if every principal subtensor has no 
nonpositive H-eigenvalue with a positive H-eigenvector. 

By Theorem 12.21 each principal subtensor of a P- or Po-tensor is also a P- or Po-tensor, 
thus has no nonpositive or negative H-eigenvalue, respectively. Then the relationship between 
P- or Po-tensors and (strictly) copositive tensors follows directly by the above lemma. 

Proposition 3.3. A symmetric P-tensor is strictly copositive, and a symmetric Po-tensor is 
copositive. 

Completely positive tensors, introduced by Qi et al. [24], are defined by sums of nonneg¬ 
ative rank-one tensors. For several vectors x& £ R" fc (k = 1,2,..., m), their outer product 

xi <g> x 2 ® • • • <8> x m is an m th -order rank-one tensor with 

(Xi 0 X 2 ® X-m) .. .im — (^l)ii (X2)i 2 ■ • ■ (^m)z m : ik — 1,2,..., 77/c, k — 1,2,..., 777. 
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When x := xi = X 2 = ■ • • = x m , we often abbreviate their outer product as x® 7 ". Then the 
(strongly) completely positive tensors are defined as follows. 

Definition 3.3 (Completely positive mm)- Let A be an m th -order n-dimensional sym¬ 
metric tensor. Then A is called completely positive (CP) if it can be written into 

r 

A = ^ u® m with Ufc £ R" for all k = 1,2,..., r; 

fc=l 

And A is called strongly completely positive (SCP) if ui, U 2 ,..., u r span the whole space R”. 

In the even-order case, it is obvious that a CP tensor must be PSD and an SCP tensor 
must be PD. Not surprisingly, this also holds that odd-order CP (SCP) tensors are symmetric 
Po- (P-)tensors. 

Proposition 3.4. A completely positive tensor is a Po-tensor, and a strongly completely 
positive tensor is a P-tensor. 

Proof. We only need to prove the case when m is odd. Let A = Y(,k=i u ?”\ where uj, £ R” 
for all k = 1, 2,..., r. Hence, A is a CP tensor. For an arbitrary x £ R”, we have 

r 

A*™- 1 = ^ Ufe (x T Ufe r- 1 . 

£=i 

Because (to — 1) is even and Ufc’s are nonnegative, the tensor-vector multiplication _4x m_1 
must also be nonnegative. Thus A is a Po-tensor. 

When Ui, U 2 ,..., u r further span the whole space, for an arbitrary nonzero x £ R", there 
must be some u*, which is not orthogonal with x, i.e., x T u^. 0. Then there is some index i 
such that Xi ^ 0 and (ufc)j > 0. Note that (ufc)j > 0 implies (*4x m-1 )j > 0. Thereby, A is a 
P-tensor. □ 

Very recently, Luo and Qi da generalized the doubly nonnegative matrices to higher- 
order doubly nonnegative tensors with all entries and all fL-eigenvalues nonnegative. Doubly 
nonnegative matrices play an important role in the relaxations of some NP-hard optimization 
problems to nonnegative semidefmite programming. Luo and Qi showed a potential similar 
effect of doubly nonnegative tensors. 

Definition 3.4 (Doubly nonnegative). A symmetric tensor is called doubly nonnegative 
(DNN) if both its entries and H- eigenvalues are all nonnegative. 

The relationship between DNN tensors and Po-tensors follows directly by Theorem 12.11 
and Proposition 13.II 

Proposition 3.5. A symmetric nonnegative Po-tensor is DNN, and an even-order DNN 
tensor is a Po-tensor. 

It is worth pointing out that in the above proposition, the assertion in the second part 
fails for odd-order tensors, as the following counterexample illustrates. 

Example 3.1. Denote A as the symmetric nonnegative 3 ld -order 3-dimensional tensor with 
am = 100, 0222 = 3, 0333 = 1, 

Oll2 = 0113 = 0122 = «133 = 1, «223 = 3, 0233 = 2.5, 

0123 = 0. 
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Note that for any nonzero x G K , (_4x 2 )i > 0 since the matrix A(l ,:,:) is positive definite. 
This implies that all H-eigenvalues of A is positive. Thus A is a DNN tensor. However, by 
setting y = (0,1, — 1) T , it follows that y|(-4y 2 )2 = —0.5 and y 2 (Ay 2 )o = — 1. Thus, A is not 
a Po-tensor by definition. 

To conclude, the three existing extensions of positive definiteness to odd orders are more 
or less involving some extra restrictions: The definition of copositive tensors constrains the 
vectors in the nonnegative orthant, and CP tensors and DNN tensors must be nonnegative 
tensor themselves. By contrast, our extension, i.e., symmetric P-tensors, admits no redundant 
conditions and best corresponds the performance of even-order PSD tensors. 

4 Structured P- and Pp-Tensors 

In this section, we will discuss the relationship between P- or Po-tensors and several important 
structured tensors. Our results show that the sets of P- and Po-tensors contain a large amount 
of structured tensors, which is another motivation for us to investigate P- and Po-tensors. 

The first kind of structured tensors is H-tensor, whose definition is based on M-tensor. 
Thus we introduce the concept of M-tensors first. M-tensor is a generalization of famous M- 
matrix and was widely studied by Zhang et al. |36] and Ding et ah [7]. Its original definition 
is as follows, where p(-) represents the spectral radius of a tensor, i.e., its largest eigenvalue 
in absolute value sense. 

Definition 4.1 (M-tensor). Let A = si — B be an m th -order n-dimensional tensor, where 
B is nonnegative. If s > p{B), then A is called an M-tensor. If s > p(B), then A is called a 
nonsingular M-tensor. 

Given an m th -order n-dimensional tensor A, the comparison tensor AA (A) of A is defined 
by 

I M( A)\ _/ l> if *1 = *2 — ' • ' = 

l i | —|o*ii 2 ...i m |) otherwise. 

The comparison tensor focuses on the absolute values of the tensor entries rather than their 
signs. H-tensor was first introduced by Ding et ah [Jj, and also studied by Kannan et ah j25] 
and Li et ah [T^]. We now introduce the definition of H-tensors. 

Definition 4.2 (H-tensor). Let A be an m th -order n-dimensional tensor. If Ai(A) is an 
M-tensor, then A is called an H-tensor. If AA( A) is a nonsingular M-tensor, then A is called 
a nonsingular H-tensor. 

As we have proven, a P- or Po-tensor owns all positive or nonnegative diagonal entries, 
respectively. Thus we focus on the H-tensors with all nonnegative diagonal entries. The 
following proposition explains the relationship of H-tensors and P-tensors. 

Proposition 4.1. A nonsingular H-tensor with all positive diagonal entries is a P-tensor, 
and an H-tensor with all nonnegative diagonal entries is a Po-tensor. 

Proof. Let A be a nonsingular H-tensor with all positive diagonal entries. Then its comparison 
tensor M.(A) is a nonsingular M-tensor. Assume that A is not a P-tensor, that is, there exists 
a nonzero vector x such that x 7 f l ~ 1 {Ax. rn ~ 1 )i < 0 for all i = 1,2 ,...,n. Hence, there is a 
diagonal matrix D with all nonnegative diagonal entries such that 

Ax™- 1 = -D-x^-b. 
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Denote V as the diagonal tensor sharing the same diagonal entries with D. Then the above 
equation is rewritten into (A+T^x” 1-1 = 0. That is, the tensor .4+2? has a zero H-eigenvalue. 

Since all the diagonal entries of A are positive, the comparison tensor of A + V is exactly 
A4(A) + V. Furthermore, because A4(A) is a nonsingular M-tensor, these is a diagonal 
matrix S with all positive diagonal entries such that A4(A)S m ~ 1 — AA^AS" 1-1 ) is strictly 
diagonally dominant. Thus the tensor M{A)S m ~ 1 + VS m ~ 1 = M((A + V)S m ~ 1 ) is also 
strictly diagonally dominant, which implies that ( A+V)S m ~ 1 is strictly diagonally dominant. 
This contradicts that A + V has a zero H-eigenvalue. Thereby, A must be a P-tensor. □ 

A tensor A is called diagonally dominant if 

| du...i | 'i! ^ [ | .. -i r n. I for all i — 1,2,... ,n, 

and called strictly diagonally dominant if 

\au...i\> [<W..i m [ for all z = 1,2,. ..,n. 

....>) 

Denote 1 = [1,1,..., 1] T as the all-one vector of size n. Then the definition of (strictly) 
diagonally dominant can be rewritten into M{A) l m_1 > (>)0. By Theorems 3 and 16 in 
(7], we know that A 1(A) is thus a (nonsingular) M-tensor, i.e., A is a (nonsingular) H-tensor. 
Therefore, we have the following corollary of Theorem ld.il 

Corollary 4.2. A strictly diagonally dominant tensor with all positive diagonal entries is a 
P-tensor, and a diagonally dominant tensor with all nonnegative diagonal entries is a Pq- 
tensor. 


Po-tensors also have some interesting examples from the spectral theory of uniform hyper¬ 
graphs (see, e.g., SUED ESI)- Actually, both the Laplacian tensor and the signless Laplacian 
tensor of a uniform hypergraph are Po-tensors. 


Definition 4.3 (Laplacian tensor [22]). Let G = (V,E) be an m-uniform hypergraph. The 
adjacency tensor of G is defined as the m th -order n-dimensional tensor A whose (zi,..., i m )th 
entry is 


au 


(m—1)! > */ {*1) ■ • • ) *rra} G E; 

0, otherwise. 


Let V be an m th -order n-dimensional diagonal tensor with its diagonal element di,,,i being 
di, the degree of vertex i, for all i £ [n]. Then C := V — A and Q := V + A are called the 
Laplacian tensor and signless Laplacian tensor of the hypergraph G, respectively. 


For even-order hypergraphs, we already know that their (signless) Laplacian tensors are 
positive semidefinite. However, we do not have this parallel results for odd-order case. Now, 
we can character this property of both even and odd-order hypergraphs via the concept of 
Po-tensors. Definition 14.31 and Corollary 14.21 indicate the following corollary. 


Corollary 4.3. The (signless) Laplacian tensor of a uniform hypergraph is a symmetric 
Po~tensor. 


Cauchy tensor is another interesting structured tensor (see HI)- Given a vector u £ 
we can construct an m th -order symmetric Cauchy tensor C as 


Ui . +u i2 -1-b Ui 


-, for all Zfe = 1,2,..., n, k=l,2,...,m, 
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where u is called the generating vector of the Cauchy tensor C. Luo and Qi Q3] proved that a 
Cauchy tensor with a positive generating vector is CP and a Cauchy tensor with a mutually 
distinct positive generating vector is SCP. Then the following corollary of Proposition 13.41 is 
apparent. 

Corollary 4.4. A Cauchy tensor with a positive generating vector is a Po-tensor. Further¬ 
more, when the entries of its generating vector are mutually distinct , it is a P-tensor. 

Recently, Song and Qi [3ll introduced B- and Bo-tensors. In matrix case, we know that 
the set of B-matrices is a subset of the set of P-matrices. Nevertheless, Yuan and You 
gave a counterexample to show that a B-tensor is not necessarily a P-tensor for the case that 
the order is even and larger than 2 under the Song-Qi’s definition. Interestingly, we can prove 
that an odd-order B-tensor is still a P-tensor under our proposed definition. We first recall 
the definition of B- and Bo-tensors. 

Definition 4.4 (B-tensor [22]). An m th -order n-dimensional tensor A = (a,i called 

a Bq- tensor if 

^ ' aa 2 ...i m fl 0 
*2 ,...,im£W 

and 

n m— 1 ^ ^C?2 J jm ) 7^ (i )...,*) • 

j2,...,i m eW 

If strict inequalities hold in the above inequalities, then A is called a B-tensor. 

We have the following theorem to construct other P-tensors from a given one in the odd- 
order case, which also implies the relationship between odd-order B-tensors and P-tensors. 

Theorem 4.5. Let A be an m th -order n-dimensional tensor, and m is odd. If A is a P- or Po- 
tensor and C is a completely positive tensor, then A+C is also a P- or Po-tensor, respectively. 
Particularly, the addition of a P- or Po-tensor and a nonnegative diagonal tensor is also a 
P- or Po-tensor, respectively. 

Proof. Denote C = X)fc =1 (ufc) m , where Ufc’s are nonnegative vectors and 


(u fc ) m = u k ® Ufc (8) • • • <g> Ufc . 

s/ ^ 

m 

Thus (A + C)x m_1 = _4x m_1 + 1 u k ( u fc x )' n 1 ■ Since to is an odd number, the second 

part Y^k=l u fc ( U J X )” ! 1 must be a nonnegative vector. From the definition of P-tensors, if A 
is an odd-order P-tensor, then there must be some index i such that Xi ^ 0 and (_4x m-1 )i > 0. 
By the above discussion, we know that ((A + C)x m_1 ), is also positive, which implies that 
A + C is also a P-tensor. 

Every diagonal tensor T> with diagonal entries di, ^ 2 ,..., d n can be decomposed into 

n 

fc=i 

where e k is the fc-th orthonormal column vector. So it is a complete positive tensor, and the 
results follow directly. □ 
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Corollary 4.6. An odd-order B-(Bo~)tensor is a P-(P$-)tensor, and a symmetric even-order 
B-(Bo~)tensor is a P-(Po~)tensor. 

Proof. For the first part, let A be an orld-order B-tensor. Song and Qi m and Yuan and 
You 02 proved that a B-tensor can always be decomposed into 

A = M+C, 

where M is a strictly diagonally dominant M-tensor and C is the summation of several 
partially all-one tensors. Since a nonsingular M-tensor is a P-tensor and a partially all-one 
tensor is completely positive, thus a B-tensor A is also a P-tensor. The second part follows 
readily from YE>\ and Proposition ld.il □ 


5 P- and P 0 -Tensor Cones 

For simplicity, we use the following notations to denote different tensor cones. The orders 
and the dimensions of the involving tensors are m and n , respectively, by default. 


Abbreviations 

Full explanations of the cones 

T 

m,n 

square tensors 

&rn,n 

symmetric tensors 

P{m , n) 

P-tensors 

Po(m,n) 

Po-tensors 

CP m ,n 

completely positive tensors 

SCP m , n 

strongly completely positive tensors 

NSPSD m , n 

positive semidefinite tensors 

NSPD mtn 

positive definite tensors 

PSD m , n 

symmetric positive semidefinite tensors 

P Bm,n 

symmetric positive definite tensors 

Zm,n 

Z-tensors 

Mm^ n 

M-tensors 

SM mtn 

nonsingular M-tensors 


The following theorem explains that Po(m, n) is the closure of P(m, n), which is the 
same relationship with CP m ^ n and SCP mtn , PSD mtn and PD mtn , M mtn and SM mtn , etc. 
Therefore, when we come to some results about P(m, n) and those results can be preserved 
by taking limits, thus those results can also holds for Po(?n,n). 

Theorem 5.1. int(Po(m,n)) = P(m,n), and thus cl(P(m,n)) = Po(m,n), where int(-) and 
cl(-) denote the interior and the closure of the corresponding set, respectively. 

Proof. On one hand, if A is a Po-tensor, then A + el is a P-tensor for all e > 0. Thus each 
Po-tensor must be the limit of a sequence of P-tensors. On the other hand, if A is the limit of 
a sequence of P-tensors {AjfcA, then we shall prove that A must be a Po-tensor. Since each 
A is a P-tensor, for an arbitrary xel" there is an index ik such that x™ 1 (Ax m 1 )i k > 0. 
Note that the series AlfcA must admit a convergent subseries (AljA Thus there is an 
index i = linq^oo i kj such that xf l ~ 1 (Ax. m ~ 1 )i > 0 with Xi ^ 0. That is, A is a Po-tensor. □ 

It is known that Po(m,n) is not convex. Counterexamples can be found in [34] for the 
case of matrices. For general high order tensors, the following proposition can also indicate 
the nonconvexity of Po(m,n). 
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Proposition 5.2. For any j £ [m] and any ij € [n], £/ie tensor <g) 
a Po-tensor; For any j £ [m] and any ij £ [n], if = 0, f/ien —e( n ) ® 

Po-tensor. 


i('m) 


JS 


j(*m) 


*s a 


Proof. For any given j = 1,2,..., m and any ij = 1,2,..., n, denote A := e^ 0 ■ ■ ■ 0 
as the rank-one m th -order n-dimensional tensor. For any nonzero x £ R ra , it follows readily 
that for any i = 1, 2,..., m, 


(4) 




-ly = / X Z lx i2 ■ • • X im > if * = *i; 

0, otherwise. 


Denote supp(x) := {i : Xi ^ 0}. 

Case 1 If *i supp(x), that is, there exists some k ^ i\ such that Xk 0. By the equation 
10]), we have x™^ 1 ^l(x m_1 ) fc = 0. 

Case 2 If i\ = supp(x), by setting k = ii, the equation 0]) implies that 

„m-1( — l \ _ / xS \ if 5ki 2 ...i m — lj 


\Ax m ~ l ) k = 


otherwise. 


Thus, for both the cases, we can always find some index k = 1,2,..., n such that Xk 0 
and x r ff~ 1 (Ax m ~ 1 )k > 0. By Definition 12.11 A is a Po-tensor. Similarly, we can show the 
remaining part. □ 

For any set C, conv(C) stands for the convex hull of C and cone(C) for the convex cone 
generated by C. As an extension of Proposition 4.1 in 03], we have the following property 
for Po-tensor cone by invoking Proposition 021 

Proposition 5.3. conv(Po(m,n))=cone(Po{m,n)) = {de T m _ n : diag(A) > 0}. 

Proof. By invoking Proposition 15.21 we have 

conv(Po(m, n)) 

= cone(Po(m, n)) 

2 cone ({E ll - lm ,ij £ [n],j £ [m]} U {-E ll ~ Zm ,ij £ [n],j G [m], dy...^ =0}) 

T a il ...i m E 11 - 11 + ^2 > 0,6^...^ G M,V*i,. ..,i m G [n] 




= \ £ + 


L ... im =0 


^ ^ Ci^imE 1 m • ^0? Cii...%m ^ 3R., Vzi,..., z m G \n\ 


:[n\ 

= {A G P mj „ : diag(A) > 0} 

On the other hand, for any Po-tensor A , we have diag(A) > 0 from Theorem 12.21 Thus, 
conv(P 0 (m, n)) C conv({A G T m ^ n : diag(A) > 0}) = {A G T m , n : diag(A) > 0}. 

This completes the proof. □ 
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The following two propositions are about the relationship between P-tensors and PSD 
tensors. 

Proposition 5.4 ([3T]). Let m be even. We have 

P ( tti , n) P — PD rnn and Po (rn, rfj P — PSD m ^ n . 

Proposition 5.5. NSPSD m , n is a convex cone in T TOi „. Moreover, 

int (NSPSD m , n ) = NSPD m , n , NSPD m , n C P(m, n), NSPSD m , n C P 0 (m, n). 

Proof. The convexity of NSPSD m}n and the assertion int (^NSPSD m ^ = NSPD m ^ n follow 
readily from Definition 13.11 And the remaining parts are obtained by invoking Proposition 

rm □ 

Corollary 5.6. P(jm,nf Pi ^ AAn,n? and Pq(tti,ti) P Z m T i — WLm.n- 

Theorem 5.7. P(m,n) P CP mtn = SCP m ^ n . 

Proof. To get the inclusion SCP m ^ n C P(m, n)PCP mj „, it suffices to show that any tensor A £ 
SCPm.n is a P-tensor. By the definition of strongly completely positive tensor, we can find 

r 

some nonnegative vectors u^, k = 1,..., r, such that A = ( u fc)™ and spanjui,..., u r } = 

k=l 

R ra . Assume on the contrary that A is not a P-tensor, that is, there exists some nonzero 
y el" such that for any i £ supp(y) := {i '■ Vi ^ 0}, 

(5) y?~\A < 0. 

When m is even, it follows readily from (|5J) that 

0> £ y i (Ay m - 1 ) i = Ay m = j2(y T ^) m >0. 

iGsupp(y) k—1 

Thus, Ay m = 0. Together with spanjui,..., u r } = R n , we have y = 0, which is a contradic¬ 
tion. When m is odd, (j5j) indicates that for all i £ supp(y), 

0 > (iy™- 1 ), = ±u™ (y T u fc ) m_1 > 0. 

k =1 


This indicates 

u[ k \y T u k ) = 0, Vi £ supp(y), Vk £ [r]. 

Thus, 

(y T u k) 2 = (y T w) = 0 ,Vk£ [r]. 

iesupp(y) 


Since spanjui,..., u r } = R", we have y = 0, which is a contradiction. Therefore, A is a 
P-tensor. The remaining inclusion follows readily from Theorem 12.11 and the definition of 
strongly completely positive tensors. □ 
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6 Tensor Complementarity Problems and Some Other 
Applications 

The tensor complementarity problem (see imuaismEa) is referred to finding some vector x 
which satisfying 

(6) x > 0, Ax" 1 - 1 + q > 0, (x, Ax" 1-1 + q) = 0, 

which is a generalization of the linear complementarity problem and a special case of the 
nonlinear complementarity problem. The P-matrices play an essential role in linear comple¬ 
mentarity problems. That is, if A is a P-matrix, then the following linear complementarity 
problem 

x > 0, Ax + q > 0, (x, Ax + q) = 0, 

has a unique solution for an arbitrary vector q. We will see shortly that if the tensor A in the 
Problem (pj) is a P-tensor, then we have some characterization of its solution set, although 
the solution is never necessarily unique. Before that, we need an important lemma from More 

01 - 

Lemma 6.1 ( [l'6j ). Let f : R" —> R n be a continuous mapping on the rectangular cone K^, 
and assume that for each x/0 in Kr, 

max Xi(fi(x) - fi( 0)) > 0. 

2 = 1,2 

If the mapping g : R” —>■ R" defined by g(x) = f(x) — f(0) is positively homogeneous, then 
for each z € R" , there is an x* >k r 0 with 

(7) f(x*) > K * z and (x*, f(x*) - z) = 0. 

Moreover, the set of x* >k r 0 which satisfy (J7]) is compact. 

Let .Kr = R” in the above lemma and denote f(x) = Ax m_1 + q, where A is a P-tensor. 
Since f(0) = q, then the definition of P-tensors and the homogeneousness of f(x) — f(0) = 
Ax™- 1 implies that f(x) satisfies the conditions of the above lemma. Therefore, we have the 
following theorem about the tensor complementarity problems with P-tensors. 

Theorem 6.2. If A is a P-tensor, then the complementarity problem ([6]) has a nonempty 
compact solution set. 

We have another class of special tensors which guarantees the feasibility of the tensor 
complementarity problem ([6]). 

Definition 6.1 (S-tensor). Let A be an m th -order n-dimensional tensor. Then there exists 
a vector x > 0 such that Ax m_1 > 0, then we call A an S-tensor. 

From Theorem 16.21 we can conclude that P-tensors are special S-tensors. 

Theorem 6.3. A P-tensor is also an S-tensor. 

Proof. Let A be a P-tensor, then from Theorem 16.21 we know that there must exist some 
x > 0 such that 

Ax" 1-1 > 1 and (x, Ax" 1-1 - 1) = 0, 
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where 1 denotes the all-one vector of the corresponding size. Then for an arbitrary e > 0, we 
have x + el > 0 and 


A(x + el)" 1 ' 1 = Ax™- 1 + O(e) > 1 + O(e). 

When e is small enough, the tensor-vector multiplication .4(x + el) m_1 > 0. This completes 
the proof. □ 

As symmetric real tensors are closely related to higher order derivatives of sufficiently 
differentiable multivariate real-valued functions, we have the following property of symmetric 
P-tensors. 


Proposition 6.4. For any given real-valued function f : 1" —>• R which is third-order con¬ 
tinuously differentiable, any given vectors £ R ra , denote V/(x) and V 2 /(x) as the gradient 
vector and the Hessian matrix of f atx respectively, and V 3 /(x) as the third-order derivative 
tensor with 


V 3 /(x)d 3 


d 3 f 


ii ,i3£M 


dxi x dxi 2 dxi 


-di 1 di 2 di ; 


If V/(x) = 0, V 2 /(x) = 0, and V 3 /(x) is a P-tensor at some point x, then for any nonzero 
d € R", we can find some sufficiently small a(d) > 0 such that f(x + a(d)d) > /(x) > 
/(x-a(d)d). 


Proof. Write the third-order Taylor expansion of / at x as follows: 

(8) 

/(x+a(d)d) = /(x) + a(d)V/(x) T d+ a ^ ) d r V 2 /(x)d+V 3 /(x)d 3 + o (a(d) 3 ||d|| 3 ) , 


By employing Corollary 13.31 V 3 /(x) is a strictly copositive tensor. This immediately implies 
that for any nonzero d £ R+, V 3 /(x)d 3 > 0. Combining with the conditions V/(x) = 0 and 
V 2 /(x) = 0, the Taylor expansion (JHJ) implies the desired assertion. □ 

Remark 6.5. By using the similar scheme of proof, the result in Proposition can be 
generalized to any odd-order case larger than 3. 


7 Conclusions 

By employing a homogenous formula, the concept of higher-order P-tensor has been intro¬ 
duced which is well-defined for all even and odd order tensors. This class of tensors has been 
shown to contain many important subclasses of tensors, including the well-known positive 
definite tensors, the nonsingular M-tensors, the nonsingular H-tensors with positive diago¬ 
nal entries, the strictly diagonally dominant tensors with positive diagonal entries, and so 
on. As an extension of PSD tensors, the Po-tensors and their relationship among other PSD 
extensions have been characterized. Additionally, many structured P- (Po-) tensors have 
been stated and many related tensor cones have been analyzed. All of these can be treated 
as refinement components for tensor theory and analysis. As an important application of 
P-tensors, the tensor complementarity problems with P-tensors have been shown to admit 
nonempty and compact solution sets. The application in multivariate real-valued functions 
has also been discussed. 
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